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Abstract 

Using a gauge-invariant formalism we derive and solve the per- 
turbed cosmological equations for the BSBM theory of varying fine 
structure 'constant'. We calculate the time evolution of inhomo- 
geneous perturbations of the fine structure constant, ^ on small 
and large scales with respect to the Hubble radius. In a radiation- 
dominated universe small inhomogeneities in a will decrease on large 
scales but on scales smaller than the Hubble radius they will undergo 
stable oscillations. In a dust-dominated universe small inhomogeneous 
perturbations in a will become constant on large scales and on small 
scales they will increase as t"^^^ , and ^ will track . If the expan- 
sion accelerates, as in the case of a A or quintessence-3ominated phase, 
inhomogeneities in a will decrease on both large and small scales. The 
amplitude of perturbations in a will be much smaller than that of mat- 
ter or radiation perturbations. We also present a numerical study of 
the non-linear evolution of spherical inhomogeneities in radiation and 
dust universes by means of a comparison between the evolution of flat 
and closed Friedmann models with time- varying q. Various limitations 
of these simple models are also discussed. 
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1 Introduction 



A third sample of observations of quasar absorption-line spectra has been 
found to be consistent with a time variation in the value of the fine struc- 
ture 'constant' between redshifts z = 0.2 — 3.7 and the present. The entire 
data set of 128 objects gives spectra consistent with a shift of Aa/a = 
[a{z) — ao]/ao = —0.57 ± 0.10 x 10~^, where ao is the present value of the 
fine structure constant |21 El II]- Extensive analysis has yet to find a se- 
lection effect that can explain the sense and magnitude of the relativistic 
line-shifts underpinning these deductions. Motivated by these observations, 
there has been considerable theoretical investigation of the cosmological con- 
sequences of varying a, jS]. Barrow, Sandvik, and Magueijo [HI [3 |H1 El UHj 
have studied in detail the cosmological consequence of an extension of the 
varying-a Maxwell equations formulated by Bekenstein |TT]. We call these 
BSBM theories. They allow us to understand the effects of the expansion of 
the universe on variations of a and to evaluate the effects of varying a on free 
fall which leads to potentially observable violations of the weak equivalence 
principle |12[ IT^ IT^ [T3] . They also allow us to investigate whether or not 
other cosmological observations are consistent with the small variations in 
a that are required to fit the quasar observations. Other theories, founded 
explicitly upon variation in the velocity of light, have also been proposed but 
their main interest is in connection with events in the very early universe 

[T7| ITHl IT^ and the problem of whether or not non-inflationary explana- 
tions for particular features of the large-scale structure of the universe are 
possible |2ni- Others have begun to consider the implications for grand uni- 
fication of coupling strengths, [21111312211221121112312311311211131, and 
astronomical probes of the constancy of the electron-proton mass ratio have 
reported possible evidence for time variation, |28j, but as yet the statistical 
significance is low. 

Barrow, Sandvik, and Magueijo [3 [3 |H1 El have shown that BSBM theo- 
ries have a number of appealing properties. They predict that there should 
be no variation of a during the radiation era in our universe and none dur- 
ing any present or late-time curvature or cosmological constant dominated 
era. During the dust era a should grow (leading to Aa/a < as observed) 
but only as log(t).This behaviour allows the quasar data to be accommodated 
without producing conflict with recent geonuclear limits on allowed variations 
of a, like the Oklo natural reactor limits of Aa/a < 10""^ [23 EOl EH E21 El 
because they are imposed at a very low effective redshift of ^ ~ 0.15, at which 
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time the universe has begun accelerating and all variations in a are damped 
out. Recent deductions of a possible upper limit of Aa/a < 3 X 10"^ at 
z = 0.45 from nuclear P decays are potentially more restrictive IHUl IH7] . 
However, it must be remembered that both these nuclear limits are derived 
from a local solar-system environment. In the absence of a theory relat- 
ing the value of, and rate of change of, a on the cosmological scales where 
quasar lines form to their values in the virialised local inhomogeneities where 
galaxies, stars and planets form, one should be wary of ignoring the possible 
corrections that must be introduced when comparing planetary and quasar 
bounds: the density of the Earth would not, for example, be a good indi- 
cator of the density of the universe. Thus, inhomogeneity is an important 
factor in the study of varying-constant cosmological theories. In this paper 
we are going to study the evolution of small, gauge- invariant perturbations to 
the exact Friedmann-Robertson- Walker solutions of the BSBM theory. The 
results of this investigation will reveal whether one must worry about fast 
growth of small initial inhomogeneities in the value of a, which would lead 
to spatial variations in the value of a that might be more significant than 
the time variations at late times. 

Bekenstein [llj generalised Maxwell's equations to include varying e and 
this theory was then generalised by Sandvik, Barrow and Magueijo jH] to 
include gravitation. In the BSBM varying-a theory, the quantities c and h 
are taken to be constant, while e varies as a function of a real scalar field ip, 
with a = e^, hence 



Our discussion is organised as follows. In section 2 of the paper we give 
the gauge-invariant perturbation equations following the development used in 
general relativity. In section 3 we specify the BSBM cosmology with varying 
a and derive the gauge invariant linear perturbation equations which couple 
the perturbations in the gravitational field to those in a and the density. In 
section 4 we solve for the time-evolution of small inhomogeneities in the fine 
structure 'constant' on large and small scales for radiation, dust, and cosmo- 
logical constant dominated expansion of the background universe. In section 
5 we extend these studies into the non-linear regime by means of numeri- 
cal solutions for fiat and closed FRW universes. The evolution of spherical 
curvature inhomogeneities in density and in a is followed by computing the 
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difference in time evolution between the FRW models of different curvature. 
These studies also reveal for the first time the behaviour of a in closed FRW 
models in the BSBM theory. These results are discussed and conclusions 
drawn in section 6. 

Units will be used in which c — H— 1; Greek indices run form to 3 and 
Latin indices only over the spatial degrees of freedom 1 to 3. The Einstein 
summation convention is assumed; a{t) is the scale factor of the background 
Friedmann-Robertson- Walker (FRW) metric and G is Newton's gravitation 
constant. 

2 The Background and Perturbations 

In what follows we shall assume that the metric of space-time deviates only 
by a small amount from a homogeneous, isotropic FRW space-time which 
is defined to be the background universe. In this case, it is convenient to 
split the metric into two parts: the background metric, and its perturbation. 
We observe that the universe is nearly homogeneous and isotropic on large 
scales (metric perturbations to the cosmic microwave background radiation 
are observed to be small ~ 10~^ << 1). 
The background line element is: 

ds^ g^^dx^dx"" = df - a^{t)^ijdx'dx> = a^{r)){dri^ - j^jdx'dx^), (3) 
where rj is the conformal time, 

dr] — a~^dt. 

We choose the background metric to be the FRW metric, so 

where k — 0, 1, —1 depending on whether the three-dimensional hypersur- 
faces of constant rj are fiat, closed or open. The Einstein equations are: 

g;; = i?;; - ^6^,r = sttGti:, (4) 

where is the Ricci tensor, R = R>^ is the Ricci scalar and is the 
total energy-momentum tensor. We shall exploit the fact that the BSBM 
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theory for varying a can be expressed as general relativity with a particular 
linear combination of energy-momentum tensors. For the moment consider 
the presence of a single energy-momentum tensor. In the next section this 
will be decomposed appropriately. 

For the background metric in equation (jSj), in conformal time, the Einstein 
equations reduce to the — equation 

= ^-^aX - ^, (5) 
where Ti. = ^ using conformal time, and the i — i equation: 

n' + n^ = ^a^T - T = t; (6) 

For the background metric the space-space part of the Ricci tensor 
Kj is proportional to 5*. Thus, for an isotropic background universe, the 
energy-momentum tensor must also be spatially diagonal, Tj oc 5*, in order 
that the Einstein equations are satisfied. Differentiating (jSj) with respect to 
T] and subtracting 2a' we get the continuity equation for matter V^T^ = : 

(iT° = -(4T° -T)c/lna. (7) 

We now introduce small perturbations around the FRW background and 
follow the gauge-invariant approach of Mukhanov [SB] • The full line element 
may be expressed as: 

ds"^ = ^'^^gf.^dx^dx'' + 6gf,^dx''dx\ (8) 

where 6g^^ describes the perturbation. The full metric has been decoupled 
into its background parts and perturbation parts: 

9t^u = ^°^9ti,y + Sgf,u- (9) 

In general, the metric perturbations can be divided into three distinct 
types: scalar, vector and tensor. Neither of the vector and tensor perturba- 
tions exhibit growing instabilities in dust and radiation universes . Vector 
perturbations decay kinematically in an expanding universe whereas tensor 
perturbations lead to gravitational waves that do not couple to the isotropic 
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energy-density and pressure inhomogeneities. However, scalar perturbations 
may lead to growing inhomogeneities which, in turn, have an important ef- 
fect on the dynamics of matter and thereby on the time and space variations 
of the fine structure constant in the BSBM theory. In the linear approx- 
imation, scalar, vector and tensor perturbations evolve independently and 
can be considered separately. In this paper we will consider only the scalar 
perturbation modes. 

The most general form of the scalar metric perturbations is constructed 
using four scalar quantities which are functions of space and time coordinates: 

^9,.- 2(x7., 

where \i represents the three-dimensional covariant derivative. 
From the above equation and eqn. Q, the line element for the back- 
ground and for the scalar metric perturbations is 



ds^ = a2(r/){(l + 2cf>)dr]^ - 2B\idx'dr] - [(1 - 2x)lij + 2E\ij]dx'dx^)}. (10) 
2.1 Gauge- invariant Variables 

Gauge-invariant variables are unchanged under all infinitesimal scalar coordi- 
nate transformations, so they are independent of the background coordinates. 
Such quantities can be constructed out of the four scalar functions $, 
and B . The simplest gauge-invariant linear combinations which span the 
space of gauge-invariant variables that can be constructed from the metric 
variables alone are: 

\(B - E')a\' a , 

a a! 
In general, a scalar quantity /(//, x) defined in the spacetime can be split 
into its background value and a perturbation /(//, x) = /o(?7) + Sf{ri,x.). 
Since, in general, 6f is not gauge invariant, we cannot use this scalar quantity 
without modification if we want to have gauge-invariant equations. Hence, 
we consider the following gauge-invariant combination: 

6f^''^ = 6f + f;,iB-E'). (12) 
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The freedom of gauge choice can be used to impose two conditions on the 
four scalar functions. The longitudinal gauge is defined by the conditions E = 
B = 0. This gauge choice has the advantage of ruhng out the comphcations 
of residual gauge modes. Also, in this gauge cj) and x coincide with the 
gauge-invariant variables $ and "if respectively. In this longitudinal gauge, 
the metric takes the form: 

ds"^ = a'^{ri)[{l + 2^)drf - (1 - 2^)-f ^^dx' dx% 

and the gauge invariants $ and \l/ become the amplitudes of the metric 
perturbations in the longitudinal coordinate system. In the case where there 
are no space-space components in the energy-momentum tensor, so Tj oc 5^-, 
we have that $ = \l/ and there remains only one free metric perturbation 
variable which is a generalisation of the Newtonian gravitational potential. 

3 Linear Theory for Cosmological Perturba- 
tions 

For small perturbations of the metric, the Einstein tensor can be written as 
GJ^ = ^^^C^ + SC^, and the energy- momentum tensor can be split in a sim- 
ilar way. The fully perturbed Einstein equations can be obtained for scalar 
perturbation modes with the line element given by ()1U|). However, these 
equations are not gauge invariant, since they contain non gauge-invariant 
quantities. In order to have gauge-invariant equations we need to replace cj) 
and X by the gauge- invariant variables $, \1/ and B — E' and to construct the 
gauge-invariant equivalents of and dT^J- we need to rewrite then as [HH] : 

5^''^ ° = <5G° + ((°)G°)'(5-E'), (13) 

5Gr^' = 5G^ + rGl~\^'^Gl){B-E% 
and analogously for STj^, 
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The components of the perturbed Einstein equations hnearised around 
small perturbations of the background are SG^^^ ^ = StcGSt!)^^^ 

5G^i = di{a^)' = ArtGa6T^^''^° (14) 
6GI = + 3n<^' + <l>{2n' + -k) = -A7TGa^6T^^'^' 

where we have already simplified the equations since the energy-momentum 
tensor that we will be considering in section 3 has no space-space components 
and so $ = \E'. 

In order to close our system of equations, we need equations of motion 
for the matter formulated in a gauge-invariant way. This requires explicit 
expressions for the energy-momentum tensor and so we must now specify 
the BSBM theory. 

3.1 The Model and the Background Equations 
3.1.1 The BSBM Theory 

The action for the universe in the BSBM theory is given by: 



S Id X\/ ~Q (^^grav ~^ ^matter ~l~ ^i/i ~l~ ^em^ ^) ; (-^^) 




where £^ = ^df^ipd^ip, a; is a coupling constant, Cem = —\f^iuf'^^^ and ip was 
defined in eqn. ([T)). The gravitational Lagrangian is the usual Cg = —j^^R, 
with R the curvature scalar, and we have defined an auxiliary gauge potential 

= eAfj^, where e{x^) = e/eo describes change in the electron charge away 
from a constant reference value Cq. The field tensor /^,^ = eF^^, = d^a^ — d^a^, 
so the covariant derivative takes the usual form, = + ie^a^. The 
dependence on e in the Lagrangian then occurs only in the kinetic term for 
e and through the F'^ = p /e^ term. 

It was shown in jj^ and that, in the context of the BSBM model, the 
homogeneous evolution of if) does not create significant metric perturbations 
at late times and the cosmological time-evolution of the expansion scale- 
factor is very well approximated by the usual power-laws found in k = 
FRW models filled with a perfect fluid. 
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Therefore, we will assume there are no major modifications in the per- 
turbed spacetime which would lead to changes of the behaviour of the per- 
turbed variables of a perturbed FRW spacetime filled with perfect fiuid. That 
is, we will assume that the energy-density perturbations and the metric po- 
tential are the same as in a FRW universe with no variation of a. Physically, 
this is to be expected for most of the evolution, although this assumption 
might break down (along with much else) on approach to initial and final cos- 
mological singularities. It is a refiection of the fact that the changes in a have 
negligible feedback into the changes in the expansion, which are governed to 
leading order by gravity. The principal effects are those of perturbations in 
the matter density and expansion rate on the evolution of a. This simplifi- 
cation will allow us to write the time and space variations of the scalar field, 
bil), as a functions of p^, 5p^, 5p^, $ and ■?/', where 5p^, bp^ and $ will 
be given by the solutions found in ref. |38j for universes with no variation of 
a; the field will be given by the solutions found previously in JU]. In order 
to find an expression for btl) in terms of these quantities we need to write the 
gauge-invariant linearly perturbed Einstein equations for the BSBM model. 



3.1.2 The Background Equations 

We vary the action p5|) with respect to the metric to obtain the generalised 
Einstein equations: 

where T^^* = ^^^^^^^isii jg w^t^ energy-momentum tensor for perfect-fiuid 
matter fields, ancf 

where = 6q is the comoving fiuid 4-velocity; ^ and T^"" are the 
energy-momentum tensors for the kinetic energy of the field ip and the elec- 
tromagnetic field respectively: 

Note, the total energy density of the electromagnetic field is the sum of the 
Coulomb energy density CPm ^'^d the radiation energy density p^ , where 
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— 1 < < 1 is the fraction of mass density of matter in the form of the 
Coulomb energy. We will then consider T^"^ ^ as a perfect fluid: 

T^l^ = {\C\Pm + Pr+Pm+ Pr)e~'^^U^U^ - (p^ + Pr) C'^^g^^ 

The propagation equation for ip comes from the variational principle as: 

[V^g'-'d^^p] = --V^e-^^C,^. (16) 

CO 

This equation determines how e, and hence a, varies with time. It is clear that 
Cem vanishes for a sea of pure radiation because Cem = {E"^ — -B^)/2 = 0. In 
order to make quantitative predictions we need to know how non-relativistic 
matter contributes to the right hand side for equation (fTH|. through Cem = 

The background equations can now be explicitly obtained: 

= SnGa^ (^p^ + {p, + \C\pJe-^'^ + '^a-^ij'^^ +a^A-3K (17) 
3n' = -AnGa^ (p„ (l + |C|e"'^) + 2p,e'^^ + 2a;a-V') + a'A 

where A is the cosmological constant. The equation of motion for the ip field 
is: 

2H^' + ^" = ^\^p^e-'^ (18) 

UJ 

The conservation equations for the matter fields, p^ and p^ respectively, 

are: 

p'^ + SHp^ = (19) 
pl + mpl = (20) 

3.2 Linear Perturbations 

The perturbed components of the total energy-momentum (T^"*'*' = T^"* + 
T^j^ + Tll^) arise from perturbations of the different matter fields which are 
time and space dependent. In particular, we have p„ p.^ + 6p^, p^. —>■ 
Pr + and ip ^ ip + 6ip. Note also that we have to perturb the fluid 
4-velocity field, so we have Ui —>■ Ui + 6ui, where i = 1, 2, 3. 
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In order to have gauge-invariant equations we need to express the per- 
turbed energy-momentum tensor in terms of the gauge-invariant energy den- 
sity, pressure, scalar field and velocity field perturbations . The gauge invari- 
ants 5p"m \ Spi^^^ and Sip^^''^ are defined in the same way as the gauge- invariant 
perturbation of a general four-scalar, see equation (fTTj) . so: 

6pl^'^ = 6p^+p'„,iB-E'), 6p^3'^ = 6p+p'{B~E'), 5^^^^'^ = 

where 5p is the perturbed pressure for a specific component and the gauge- 
invariant three- velocity is given by |38j : 

5uf^ = 5ui + a{B - E%. 

The physical meaning of the quantities which enter the gauge-invariant equa- 
tions is very simple: they coincide with the corresponding perturbations in 
the longitudinal gauge. From now on we will drop the superscript {gi) since 
we will always be dealing only with gauge-invariant quantities. 

We can now write the gauge- invariant linearly perturbed energy-momentum 
tensor: 



^yma^erO ^ ^^^^ (2I) 

^j^rnatterO ^ {p^ + p^) a-^Ui, (22) 

^rj.matteri ^ -5p„,5], (23) 

dTt'' = ua-^[i,'5^' -^'"^y (24) 

dTf" = ua-^^'di5^, (25) 

6Tf' = uja-^(i;''^-i;'6i)'^6), (26) 

5Tr' = (|C|5p^ + 5p.)e-^'^-2e-2MC|p„ + p.)5V', (27) 

5Tr' = e-^H\C\p^ + p,+Pm+Pr)a~Hu,, (28) 

= (25^(p„+p,)e-2^-(5p„ + (5p,)e-2'/')5}. (29) 



We have assumed there are no anisotropic stresses in the energy-momentum 
tensors and we have considered only adiabatic perturbations; that is, we 
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consider the pressure perturbations to depend only on the energy-density 
perturbations. In the dust and radiation cases this means Pm = 0, Spm = 0, 
or pr = and 6pr = ^Spj., respectively. 

The fully perturbed gauge-invariant Einstein equations can be obtained 
from (|THl using the expressions above for 6TI^ f|2HI27|) . We have 

$ (^-37^2 + 3k + AGmuilj'^^ + 7^$ - AGnujip'Silj' - 3W = 

AGna'e-'^ [{e'"^ + |C|) 6p„, + 6p,. - 25^ (p, + |C|p^)] (30) 



7i;v2$+v'$' = AG7Tu^'y^6ilj-^ae-^^ [4p, + 3 {e^^ + |C|) p^] V5m (31) 



$ (n"^ + 2n' -k + agttuj-^'^^ + sn^' + $" = 

^ [a^e-^^ {6p, - 2pM) + 3cuV^'(5V^'] (32) 

It is useful to write the perturbed energy-momentum conservation equa- 
tions for each component: 



2\C\ 

Sip" = ^e-^^a^ [5p„ + 2p„ ($ - dip)] - 2nSij' + W^Sip + 4^'$', (33) 

5pl = ^6p, {en - 3ij') + p, (2vSu - 3Sij' - 6$')] , (34) 

Sp^ = -SnSp^ - p„ (v5u - 3$') . (35) 

From these expressions it is clear that perturbations in a are sourced by 
perturbations in the dust, but not vice versa. Hence we expect that, in a dust- 
dominated universe with varying a, the cold dark matter perturbations will 
behave as in a dust-dominated universe with no varying a. Notice however 
that the same cannot be concluded so easily for a radiation-dominated era, 
since there is a source term in equation proportional to Sip' . We expect 
this term to be negligible at large scales, but that might not be the case on 
small scales. 
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The gauge-invariant perturbation for a is given by eqn. as 



a 

From equations ((201), (ED), (1221) , using (fT7|) to simplify, we obtain the general 
form for the perturbation to the scalar field which drives variations in 
the fine structure 'constant', as 

(5p„ - 2p„$)] + 36^^ 2$ (m^ -k- AGiroot/j'^^ - 7^$ + 6K$' + $" 

This is a gauge-invariant expression for 6tp, written as a function of the 
gauge-invariant quantities $, Sp^, and 5p^. 



4 Evolution of the Perturbations 

In a spatially flat universe, k = 0, the general gauge-invariant expression for 
Sip becomes: 



(5p^ - 2p^$)] + Se"^^ $ (^6H^ - SCttw?/^'^) - + 6K$' + $" ) . 

It was shown in |U] and ^0] that the cosmological evolution of the metric 
scale factor is unchanged (up to very small logarithmic corrections) to leading 
order by the time evolution of ip. The dominant effect is that of the evolution 
of the scale factor on the evolution of through its propagation equation. 

It is known that perturbations of massless scalar fields, or scalar fields 
with a very small mass are negligible with respect to perturbations in the 
matter fields and the gravitational potential. Guided by this, in order to 
obtain the evolutionary behaviour for — , we will assume that the matter field 
perturbations, 6p^ and Sp^, and the metric perturbations, $, are unaffected 
by the ip perturbations to leading order. We will assume that these three 
quantities will therefore be the same to this order as they are in a fiat FRW 
universe filled with barotropic matter and a minimally coupled scalar field. 
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These assumptions are valid if the energy density of ip is much smaller than 
the energy density of the matter fields, so $ will be driven only by the matter 
perturbations. If we examine the perturbations in the non-linear regime it is 
confirmed that ^ -C ^(see figure (jU). 

4.1 Radiation-Dominated Universes 

In any radiation-dominated era in which the expansion of the universe is 
dominated by relativistic particles with an equation of state p = |p, we can 
neglect the non-relativistic stresses in the universe, in particular, the cold 
dark matter and the cosmological constant since » » pj^ to a good 
approximation. If we assume that = = = k and 5p„ = 0, the 
background equations of motion give the usual conformal time evolution for 
the scale factor and the energy density of the radiation: 

Pr = Pro^^^e^'^ a = ^/87iGp,^^r] 

where p^^ is a constant. 

The perturbations in the barotropic matter fluid and the potential $ 
come from the equations (jHUI) . ^T\ . The usual flat FRW solutions with 
constant a are obtained by setting the terms proportional to ip to zero, 
fl^ . The general solution of these equations can be obtained by expanding 
the physical quantities in terms of the eigenfunctions of the operator 
(where —k^ denotes the eigenvalue of this operator) and solving for each 
mode separately. Resuming the terms, the general solutions of the linearly 
perturbed equations for the potential $ and the barotropic matter are: 



^ = V '^{[wi]cos{wr]) — sm(wr])]Ci + [wr]sm{wr]) + cos{wr])]C2}e^^^ (38) 



and 



Sp^ 

Pr 



Ci{riw 



1 - ^ {r]wf 



cos^rjw) + [{rjw)'^ — l] sm{riw)} (39) 



+{ [l — {riw)'^] cos{riw) + rjw 



1 - ^ {T]wf 



sin(r/u;)}C2)e'''", 



where we have expanded the general solution in plane waves since we are 
assuming a spatially flat universe; Ci and C2 are arbitrary functions of the 
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spatial coordinates; k is the wave vector mode and w = k/\/3. Note that 
these quantities are all expressed in their gauge- invariant format. Finally, to 
calculate the explicit time dependence of Sip, we need to use the background 
solution for ip: 

^ = ^log(8iV) + ilog(^r/2)_ ^4Q) 

This was found in fS^ and ^01 for a radiation-dominated universe, where 
since C < in the magnetic energy dominated theories 
considered by BSBM. It is important to notice that in universes with an 
entropy to baryon ratio {S ~ 10^) like our own, ip does not experience any 
growth in time |7j. The constant term on the right-hand side of fj40|) domi- 
nates the solution for ijj{r)) throughout the radiation era. Numerical solutions 
confirm this freezing in of '0, and hence of a, during the radiation era, [Hj. 



4.1.1 Large-scale perturbations in a radiation-dominated era 

In the long-wavelength limit {wrj << 1) where the scale of the perturbation 
exceeds the Hubble radius, we can neglect spatial gradients. So, in this limit, 
becomes: 



a 



2 a 21] 

From this expression we can see at that on large scales the inhomogeneous 
perturbations in a will decrease as a power-law in time. This behaviour agrees 
with the one found in [Oj by other methods. 



4.1.2 Small-scale perturbations in a radiation-dominated era 

On scales smaller than the Hubble radius [wr] » 1) the dominant terms 
are proportional to V^Ci and V^C2, so the asymptotic behaviour for — will 
be: 



X Sex X 

Sip = — ^ oc — w rcos(?7w)V^Ci + sin(?7w)V^C2l e'^"" + (41) 
2 a 2 

+ — rsin(?7u;)V^C7i - cos(r7w) V^Cal e^^"" 
2rj '- 

On small scales we can see the perturbations on a will be oscillatory. This 
behaviour is new and does not coincide with the ones found in [Hj. 
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4.2 Dust-Dominated Universes 



In the case of a flat dust-dominated universe, filled with a p = fluid, we can 
assume = = = k and 5p^ = 0. Again, in order to obtain an exphcit 
expression of 6ip from equation (jHTj) . we will assume that the matter-field 
perturbations, p^ and 5p^, and the metric perturbation, $, are not affected 
by the ip perturbations to leading order. Thus, we will assume that these 
functions behave as in a perturbed flat FRW dust universe. 

In general, for a flat dust universe we have the following background 
solutions: 

_ -3 _ ^TrGp^p 2 

Pm ~ PniQ^ -I ^ ~ 2 V -I 

where p^^^ is constant. 

As before, we can calculate the most general gauge-invariant solutions 
for the energy-momentum perturbations 5p^, and for the potential $, as- 
suming that the ip field does not affect them to leading order, so their time 
dependences are p?H] : 

$ = Ci + Csr]-^ (42) 

and 

^ = ^ [W^'C, - 12C,) + [r^'y'C, + I8C2) r'] , (43) 

Pm ^ 

where Ci and C2 are arbitrary functions of the spatial coordinates. 

Once again, we need the background solution for ip, in order to calculate 
dip as an explicit function of the conformal time. We use the asymptotic 
solution 

^ = llog[2Ariog(^r/3)], 

which was found in , as an asymptotic approximation for a dust-dominated 
universe which is in good agreement with numerical solutions, and where, as 
above, = — — na^ > is a constant. 

' uj im 

4.2.1 Large-scale perturbations in a dust-dominated era 

On scales larger than the Hubble radius we can neglect the terms propor- 
tional to 9jCi,V^Ci, V^C2 and diC2] so in this limit we have the following 
asymptotic behaviour for the non-decaying mode: 
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= l'^ « -2C. - (44) 
2 a 111(77) 

Therefore, on large scales the inhomogeneous perturbations of a will not grow 
in time by gravitational instability. This behaviour can be understood with 
reference to the general evolution equation for ip in Friedmann universes. 
On large scales, where spatial gradients in ip can be neglected with respect 
to its time derivatives, we may view inhomogeneities in density and in ijj 
as if they are separate Friedmann universes of non-zero curvature {k ^ 0). 
The growth of inhomogeneity can be deduced by comparing the evolution of 
a in the k ^ universes with those in the k = model (a more detailed 
numerical study of this model will be presented in section 5 below). In effect, 
this uses the Birkhoff-Newton property of gravitational fields with spherical 
symmetry. We note that the (jl8j) evolution equation has the simple property 
that cannot have a maximum because ip > when ip > because > 
, 0. This result holds irrespective of the value of k = 0. Thus ip and a will 
continue their slow increase in both over-densities, under-densities and the 
flat background until we reach scales small enough for spatial derivative to 
come significantly into play. This has the important consequence that we do 
not expect large spatial inhomogeneities in a to have developed. However, it 
should be noted that the sensitivity of the observations of varying-a effects in 
quasar spectra is sufficient to discern variations in redshift space smaller than 
O(10~^), which is of the same order as the amplitude of density fluctuation 
on very large scales in the universe. 



4.2.2 Small-scale perturbations in a dust-dominated era 

On scales smaller than the Hubble radius the dominant terms are those 
proportional to V^Ci and V^C2, so the asymptotic behaviour for the growing 
mode will be: 

5V oc ^^'v^Ci (45) 

This also shows that perturbations of a will grow on small scales. This re- 
sult is a product of the assumption that on small scales the universe can be 
considered as being filled by an homogeneous and isotropic fiuid, however we 
know this is not true below the scale where gravitational clustering becomes 
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non-linear. On these small scales we also have to worry about new conse- 
quences of inhomogeneity which have not been included in our analysis. For 
example, the constant parameter N = —"^p^a^ oc (ilm will vary in space 
due to inhomogeneity in the background matter density parameter and in 
the dark matter parameter — 1 < C ^ 1- We have assumed that C < for the 
cold dark matter on large scales in order for the cosmological consequences 
of time-varying a to be a small perturbation to the standard cosmological 
dynamics. But on small scales the dark matter will be baryonic in nature 
and so C > there. Hence, we expect ( to be significantly scale dependent 
as we go to small scales. This behaviour will be investigated elsewhere along 
with the problem of the clustering of inhomogeneities in p and a. 



4.3 Accelerated Expansion 

In an era of accelerated expansion, a > 0, as would arise during inflation or 
during a A- or quintessence-dominated epoch at late times, we can consider 
the scale factor to evolve as a power-law of the conformal time as a = r/"", 
where tt, > 1 and rj runs from —oo to 0. The case of a = ?7~^ corresponds to 
a A-dominated epoch. 

As in the previous sections, we will assume that all the other matter 
components which fill the universe will behave exactly as in a perturbed 
FRW universe with no variations in a. We also assume that neither of the 
dust and radiation perturbations will affect the behaviour and evolution of 
Sijj, and that these perturbations are negligible with respect to the A stress 
driving the expansion, so we will consider 6p,^ = and 6p^ = 0. 

It was found in [Bj and lOj that in universe which is undergoing acceler- 
ated expansion, the asymptotic solution for ip is a constant, so we will assume 
that ip = ip^, in the background, where ip^ is a constant. Thus, equations 
and (jnS) become: 

8Gn\C\5^Pp^ , ^ 3n (r/$ - n<l>') _ ^ 



n (2 + n) $ + (^r]^" - 3?2$' j = 

where we have also considered = 0, but p^ ^ because of the coupling 
with ip in the equation of motion of the scalar field. Note that if we had also 
set = here, we would have imposed a no a-variation condition: 6ip = 0. 
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Integrating the last equation, we obtain 



where Ci and C2 are arbitrary functions of the spatial coordinates. Note 
that since n > 1 in accelerating universes we see that $ will decay in time 
as 77 — > 0. Prom this solution for we obtain: 

5^1J = '-^ = -^ n 3n (l + n- Vl + n(5n-2))Ci+ (46) 

3n(l + n+ Vl + n(5n-2)) r/V^^+^^^^Cs + 2r)^S7^C, + 2r72+Vi+"(5n-2)^2^^j 

which is a decaying function of the conformal time when n > 1 and a constant 
when n = 1. Thus, in accord with the expectations of the cosmic no hair 
theorem, the universe approaches the FRW model and a is asymptotically 
constant at late times. 



4.3.1 Large-scale perturbations during accelerated expansion 

On scales larger than the Hubble radius we can neglect the terms proportional 
to the spatial derivatives; so in this limit we have the following asymptotic 
behaviour: 



5ip oc 



3g2V'oo72 n-3-Vl+7i(5n-2) 
-7] 2 



16Gp^o7r|C| 
+ (l + n+ Vl + n(5n-2)) ^V'l^M^C', 



(l + n- ^yl + n{5n-2)^ Ci 



(47) 



Therefore, as expected, on large scales during an accelerated era inhomo- 
geneities in a will decrease on time when n > 1 and will be a constant when 
n = 1. 



4.3.2 Small-scale perturbations during accelerated expansion 

On scales smaller than the Hubble radius the dominant terms are the ones 
proportional to V^Ci and V^C2, so the asymptotic behaviour will be 
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5ip oc 



V2C,+^Vl+«(5«-2)v2C2 



(48) 



This confirms that perturbations of a, as ?7 — > 0, will decrease on small scales 
when n > 1 and will be constant when n — 1. 



4.4 Summary of behaviour 



In Table 1 we summarise the time evolution of small inhomogeneities in a 
found under different conditions in this section. 



Universal equation of state 


Time Evolution of the perturbations 5i/j 


Large scales 


Small scales 


Radiation- dominated epoch 
p—\p, a OCT] 


Decaying 


Oscillatory 


Dust- dominated epoch 

p = 0, a cc rf 


Constant 


Growing 


Accelerated expn a cc ri^'^,n > 1 
A-dominated 
p — —p, a oc Tj"^ 
Power-law acceleration, n > 1 

p — wp, w < 0, a oc ?7~" 


Constant 


Constant 


Decaying 


Decaying 



Table 1: Time evolution of small inhomogeneities in a. 



5 The Non-Linear Regime 

In order to study the evolution of inhomogeneities in a beyond the domain 
of linear perturbation theory we need to use a different model. The simplest 
approach is to confine attention to spherically symmetric inhomogeneities. 
This will be done by comparing the solution of the BSBM theory for a in a 
closed (k — 1) universe, with the solution for a in a fiat {k, = 0) universe. 
We are assuming a Birkhoff property for the BSBM theory so that we can 
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Figure 1: The evolution of a{i]) and a{i]) for radiation- dominated universes with 
K = (dashed) and k = 1 (solid). 




O.E 1 1.5 2 2.5 3 

11 



Figure 2: The evolution of -^/-^ vs rj in radiation- dominated universes. 

treat the perturbation as an independent closed universe. This is a standard 
technique in general relativity which was first used by Lemaitre 

We define the alpha 'over-density perturbation' (which is not necessarily 
small) by 

a 0^=0 

where is the solution of equation (fTBj) for a universe with curvature k. 
5.1 Radiation- dominated Era 

The scale factor for a radiation-dominated closed [k = 1) FRW universe is 
given by a = sin(r7); for a fiat [n = 0) FRW universe the normalised scale 
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Figure 3: The evolution of a[r}) and a{r]) for dust- dominated universes with k = 
(dashed) and k = 1 (solid). 

factor is given by a = 77. 

The evolution of a can be seen in FigureQalong with the evolution of a for 
K = 1, 0. As expected, we can see there is no difference in the evolution of a 
at early times and a oc as it was found in ilOj. When the difference between 
the scale factors of the two universes becomes significant, the behaviour of 
aK=i begins to deviate from that of aK=o- We see that a^=i clearly grows 
faster than 0^=0 . The difference in the growth rates become very significant 
near the expansion maximum of the bound region [r] ^ tt). However, after 
this time our assumption that the background is not affected by changes of ijj 
in the cosmological equations that describe the background universe breaks 
down, since the kinetic energy of the scalar field will diverge and can no 
longer be neglected in the Friedmann equation. We expect the behaviour 
near the final singularity to be similar to the kinetic-dominated evolution 
near the initial singularity discussed in ref. jTj. 

From Figure El we see that the variations in a will become increasingly 
important as f] approaches the second singularity. Notice that although there 
is a considerable growth in the perturbations in a, when we compare them 
with the perturbations in the radiation, they are not as significant as can be 
observed from the evolution of the ratio versus 77 in Figure |21 

5.2 Dust-dominated Era 

During the dust-dominated phase of a closed universe (k = 1) the normalised 
scale factor is given by a = 2 (1 — cos^r])), while for a flat universe {k = 0) 
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10^ 10"^ 10"' 10° 10°' 10°^ 10°° 10°' 10°° 10°° 

log Ti log 11 

Figure 4: The evolution o/ log — (solid) and log (dashed) vs log(7/) for a 
dust- dominated universe. The dotted lines correspond to evolution in the linear 
perturbation solution. 



the normalised scale factor is given by a = rj'^. Integrating equation ()16|) for 
both cases we obtain the evolution of a for both cases. 

The evolution of a{ri) can be seen in Figure |21 along with the evolution 
of a{ri) for k = 1,0. As in the radiation case, we can see there is little 
difference in the evolution of a at early times, since the scale factor for the 
closed model evolves very similarly to the flat one for 77 << 1 and a oc In(ry). 
The differences between aK=i and aK=o cases start to appear when rj ^ 
1, when the nonlinear regime commences. These differences become more 
accentuated near the second singularity, but once again this is the region 
where our approximations break down. 

Notice that although there is a considerable growth in the perturbations 
in a, when we compare them with the perturbations in the cold dark matter, 
they are not as significant as can be observed from the evolution of log(^) 
and log(^^) vs. log(?7) in Figure |3] In this case perturbations in a are even 

Pm 

less significant than in the radiation case. We note that the fact that the 
linear regime is a very good approximation at early times, since as can be 
seen from Figure EJ — tracks oc t^^^ at early times. From the detail of 
Figure |4| we see that variations in the cold dark matter will start to occur 
before than variations in a, and is at least three orders of magnitude 

' Pm ^ 

bigger than — . 
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6 Conclusions 



By applying the gauge-invariant formalism of ref. j^j to the simple BSBM 
theory of varying a we have determined the evolution of small inhomo- 
geneities in baja in the presence of small adiabatic density inhomogeneities. 
To leading order, the evolution of the perturbations to the expansion dy- 
namics and the matter and radiation content of the universe behave as in 
cosmological models with constant a and we can to determine the behaviour 
of small inhomogeneities in baja in the gravitational fields created by the 
density and metric perturbations (see also refs. P and [TU]). 

In a flat radiation-dominated universe we find that inhomogeneous per- 
turbations in a will decrease on large scales while on scales smaller than the 
Hubble radius they will undergo bounded oscillations. In reality, we expect 
dissipation of the adiabatic fluctuations to occur by Silk damping and small- 
scale fluctuations in a will also undergo decay as their driving terms damp 
out. However, while the exact solution for the evolution of a is a linear sum 
of a constant and a slow power-law growth the power-law evolution does not 
become dominant by the end of the radiation era in universes like our own 
with entropy per baryon O(IO^). 

In a flat dust-dominated universe small inhomogeneities in a will become 
constant on large scales at late times while on small scales they will increase 
as t^l'^ . In reality, the small-scale evolution will be made more complicated 
by the breakdown of the assumptions underlying the perturbation analysis 
and the development of local deviations from the FRW behaviour. 

In an accelerated phase of our universe, as is the case for an early infla- 
tionary epoch, or during a A- or quintessence-dominated late phase of evo- 
lution, we show that inhomogeneous perturbations in a will decrease on all 
scales. This result complements the earlier discovery 0, 05 that a tends to a 
constant with exponential rapidity in Friedmann universes that become dom- 
inated by a cosmological vacuum stress. Any pre-existing inhomogeneities 
will be frozen in but their scale will be exponentially increased by the de 
Sitter expansion. 

These perturbative results are quite good approximations when we con- 
sider large scales, but are expected to break down when extended to small 
scales, where non-linear effects come into play and local deviations from 
isotropy and homogeneity of the matter content are significant. We note 
that the background solutions for about which we have linearised the per- 
turbations of the Einstein equations, are solutions which describe the time 
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evolution of on a 'standard' FRW background. Our neglect of the back- 
reaction of the perturbations on the background expansion dynamics is a 
good approximation up to logarithmic corrections. 

When we examined the non-linear evolution of spherical inhomogeneities 
by means of a comparative numerical study of flat and closed Friedmann 
models we found that perturbations in a remain almost negligible with re- 
spect to perturbations in the fluid that dominates the energy density of the 
universe at the same epoch. Comparing the fiat with the closed solution 
for a, we concluded that in both the cases of a radiation or dust-dominated 
epoch, a will 'feel' the change of behaviour in the scale factor, and the per- 
turbations — will grow in time. The early-time behaviour of the non-linear 



changes in proportion to at early times. 

We have provided a detailed analysis of the behaviour of inhomogeneous 
perturbations in a and its time variation on large scales under the assumption 
that the defining constant of the BSBM theory, is constant and negative 
in sign. In reality this assumption will break down on small scales. The 
negativity of the effective value of C, requires that the cold dark matter is 
dominated by the magnetic rather than the electric field energy (see also the 
discussion of ref and by Bekenstein jlHI)- However, on sufficiently small 
scales the dark matter will become dominated by baryons and the sign of the 
effective C will have to change sign. Overall, there will also be a gradient in 
the value of |C| reflecting the scale dependence of the relative contribution 
of dark matter to the total density of the universe. We have not included 
these effects in the present analysis. They would need to be included in 
any detailed analysis of the small-scale behaviour of inhomogeneities in a. 
This is an important challenge for future work because it would enable the 
quasar data on varying a to be compared directly with the limits from the 
Oklo natural reactor |3HI22] and Rhenium- Osmium abundances in meteorites 
EZl . At present the relation between the cosmological and geonuclear 
evidences is unclear because the latter are derived from physical processes 
occurring within the cosmo logically non-evolving solar system environment. 

Variations in a also affect the cosmic microwave background radiation 
spectrum and anisotropy in different ways, but the effects must be disentan- 
gled from allowed changes in other cosmological parameters that can con- 
tribute similar effects. These changes in the microwave background with a 
left as a free constant parameter were analysed in [41J using the new WMAP 
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solutions confirms the linear behaviour found in section 4. In 



a 
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data. They are far less sensitive that the many-muhiplet analyses of quasars 
aX z = 0.5—3 [1112113111, although they derive from higher redshifts, z < 1100. 
These studies can accommodate constant and varying a but up to a level 
that would be too large to be consistent with the quasar data and the slow 
time-evolution of the theory with time- varying a described in this paper. 
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